A stochastical model for periodic domain structuring in ferroelectric crystals 
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A stochastical description is applied in order to understand how ferroelectric structures can be 
formed. The predictions are compared with experimental data of the so-called electrical fixing: 
Domains are patterned in photorefractive lithium niobate crystals by the combination of light- 
induced space-charge fields with externally applied electrical fields. In terms of our stochastical 
model the probability for domain nucleation is modulated according to the sum of external and 
internal fields. The model describes the shape of the domain pattern as well as the effective degree 
of modulation. 
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The importance of ferroelectric optical materials is still 
contrasted by a lack of theoretical understanding. Com- 
paring textbook knowledge with real experimental results 
yields many discrepancies. Some examples for the promi- 
nent ferroelectric materials lithium niobate and lithium 
tantalate: Domain walls can be pinned and can bow [l|, 
the coercive field can be reduced by an order of magni- 
tude by small changes in the lithium content @ , and the 
materials exhibit a memory effect for recent domain in- 
versions [|[ . The reason that it is so hard to predict and 
to model ferroelectric domain inversion and patterning 
results from the influence of defects (vacancies, impuri- 
ties, ions on wrong lattice sites, etc.). The systems are 
too large to make first-principle calculations that con- 
sider these effects, although first attemps already yield 
impressive fundamental predictions (4j. 

The simplest periodic structure, that is ideally suited 
to compare theoretical and experimental data, is an el- 
ementary grating. Such periodically poled crystals are 
also of relevance for a variety of applications: Frequency 
conversion [5(, parametric oscillation [6(, the generation 
of terahertz radiation [3] as well as high-speed electro- 
optical switching Conventional periodic poling is 
achieved by applying high electrical fields with structured 
electrodes However, quality and size of domain struc- 
tures fabricated by this method seem to be limited. The 
reasons are not clear because of insufficient theoretical 
modelling 

One possibility to describe ferroelectric domain rever- 
sal is the theory of Avrami, Kolmogorov, Johnson, and 
Mehl on first-order phase transitions adapted by Seki- 
moto to a form usable for our purposes [HI- This 
theory relys on a purely stochastical description. Ac- 
cording to it, domains nucleate and grow independently 
from each other, and within its framework the probabil- 
ity for domain inversion can be calculated. Usually the 
domain switching is studied by obtaining the time depen- 
dence of the inverted volume fraction during the poling 
process. It is an open and demanding question whether 



this theory can be applied to domain patter ning as well. 

In the present work electrical fixing [l2], [l3[ [3, Ell 
is used to compare stochastical modelling of domain pat- 
terning to experimental data: Space-charge fields are cre- 
ated in a crystal via the photorefractive effect and 
additionally an external electrical field is applied. Ac- 
cording to the stochastical model a spatially modulated 
electrical field leads to a modulated probability for do- 
main inversion. The dependence of the domain grating 
quality on the underlying probability density can be cal- 
culated by adapting the theory of Avrami, Kolmogorov, 
Johnson, and Mehl to the case of a spatially sinuosoidally 
modulated probability for domain nucleation. Sekimoto 
gives the following general formula for the probability 
w that a point r has an inverted spontaneous polariza- 
tion: 



dt' d 3 r'I(r',t') 1 - D(f, t; P, t') 



w{r, t) = 1— exp 



Here I(r',t') denotes the probability density for domain 
nucleation at position r' and time t' and V denotes the 
whole crystal volume. The function D(r, t; r', t') is zero if 
such a domain contains r at time t and is one otherwise. 
The time dependence in D and / is introduced to model 
the temporal evolution as domain nucleation probability 
and domain sizes may vary with time. As we are not 
interested in the temporal evolution here, we replace the 
time dependence by the dependence on domain length I . 
Thus it follows: 

w(r) = 1-exp {- j dl' J d s r'l(P,l') l-D(r,P,l') 

(1) 

Here I(r', I') is the probability that a domain located at 
r' with length I' can be found. 

We now assume that due to a periodic space-charge 
grating the probability for domain nucleation is periodi- 
cally modulated. Thus we write: 



I(P,l') = [a + f3sm(Kr' z )]nt{l'). 



(2) 
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Here ni(l') describes the size distribution of domains. 
So (3/a is a measure for the strength of the influence 
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of the space-charge grating on domain nucleation. Here 
K = 2n/A is the spatial frequency of the space-charge 
grating. The domain shape is assumed to be cylindrical 
with length 21 and area A(l) and to be centered around 
r' which is described by the function D . Instead of a 
modulated probability for domain nucleation, a modu- 
lated domain-wall growth can be assumed by adapting 
I(r' , V) . However, within the framework of our methods 
both assumptions lead to very similar predicitions, hence 
the latter is not investigated within this article. Insert- 
ing equation (f2]) into equation {T]) and using the specific 
domain form, we get: 

w(r) = 1-exp ( - 2a dlrii(l)A(l)fj 

x exp f - 2sm{Kr z )— [ dlm(l)A(l)sm(Kl) 
v K Jo 

We define 7 = 2a Jl°° dlni(l)A(l)l and ij = 
2(J3/K) J °° dlni(l)A(l)sm{Kl). The former is a mea- 
sure for the volume of inverted spontaneous polarization 
whereas n = rj/j measures the modulation degree of the 
domain grating (0 < k < 1). Expanding the second ex- 
ponential function in the last expression for w up to the 
first order we get as an approximate expression for w. 

w(f) = 1 — [1 — 7K sin (Kr z )\ exp (—7) . 

The area fraction q with inverted spontaneous polariza- 
tion (degree of poling) can be found by taking the spatial 
average of w(r): 

1 f A 

q= Aj dz l( z ) = 1 ~ ex ^(~V ■ 

The degree of modulation of the domain grating Aq can 
be calculated as the first Fourier coefficient for K: 

Aq = jk cxp(— 7) = — k(1 — q) ln(l — q) . (3) 

Equation [3] expresses the domain grating modulation 
solely as a function of the degree of poling. For no or com- 
plete reversal of the spontaneous polarization no grating 
at all can be found. A maximum is found in between. Nu- 
merically evaluating equation [3] for several assumptions 
yields that for the height of the curve maxAg w 0.4k 
roughly holds. 

To test these predictions experimentally, electrical fix- 
ing experiments are performed. Lithium niobate crystals 
doped with 0.05 mol% iron from Deltronic Crystal Indus- 
tries are used. To reduce the coercive field a vapor trans- 
port equilibration treatment is applied [l7| . The crystals 
are put into a special mount which is placed in a glass bin 
filled with silicon oil to prevent electrical breakdowns. 

The optical setup consists of a detuned Mach-Zehnder 
interferometer, which allows to illuminate the crystal 
with a periodic interference pattern with grating peri- 
ods of 5 to 100 (im. Light of the wavelength 488 nm from 



an Argon-ion laser is used. The illumination creates a 
space-charge field which induces an index-of-refraction 
grating via the electro-optic effect. This photorefractive 
volume grating is detected by Bragg-diffraction with light 
of the wavelength 633 nm from a HeNe-laser. From the 
amount of diffracted light the amplitude of the modu- 
lation of the index-of-refraction and hence of the space- 
charge field can be calculated [18]. A detailed description 
of the setup can be found in reference [191 ]. 

The experiments are performed as follows: A space- 
charge field is written by illuminating the crystal with 
an interference grating. Either during or after the record- 
ing a voltage is applied to the crystal along the z-axis. 
The polarity of the voltage is chosen such that the elec- 
trical field supports domain inversion. The strength of 
the space-charge field is monitored during the whole ex- 
periment. When the space-charge field reaches steady 
state, the illumination is stopped and the external elec- 
trical field is shut off. After homogeneous illumination 
of the crystal no diffracted light is detected. Next the 
crystal is repoled to its single-domain state by applying 
an electrical field without any illumination. The degree 
of poling is determined by integrating the current during 
the poling process. This value is normalized by twice the 
spontaneous polarization Pg. 

The revealed index-of-refraction grating provides in- 
formation about the achieved degree of modulation of 
the domain walls. It is studied for different strengths 
and application times of the external field during the 
recording process. In Fig.Q] the revealed space-charge 
grating strength is plotted against the degree of poling. 
In Fig.[T](a) a space-charge field is written in the crys- 
tal before the external field is applied. The experimental 
parameter, which is varied, is the time this field is ap- 
plied. In Fig.[T](b) the space charge grating is written 
in the presence of an externally applied field. Here the 
recording time is fixed but in order to change the degree 
of poling the external field is varied. 

As can be seen, the strength of the revealed field be- 
comes zero for q — > as well as for q — > 1 and reaches a 
maximum in between. In Fig[T](a) the dependence of the 
strength of the reappearing space-charge field is shifted 
to higher g-values compared to the results obtained with 
varying fields in Fig.[T](b). 

Comparing theoretical predicitions with experimental 
data for the strength of the revealed field under dif- 
ferent fixing conditions, qualitative agreement can be 
found. The dependence on the degree of poling q with 
different writing times as in Fig.[T](a) is very close to a 
— (1 — x)ln(l — x) dependence as predicted by Eq. 
The shift to lower g-values, when the field is varied, is 
understandable: Different external fields do not only in- 
fluence q but also change the relative influence of the 
space-charge field and thus (3 /a. 

Experimental data for different grating periods is 
shown in Fig.O For period lengths smaller than 15 /im 
the strength of the reappearing space-charge field quickly 
drops to zero. However, for large period lengths exceed- 
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FIG. 1: The index-of-refraction modulation An is plotted ver- 
sus the degree of poling q. The period A of the fixed grating 
is A = 15.7pim. In Fig.Q](a) the degree of poling is varied by 
applying an external electrical field for different times. The 
solid line is a plot of Eq. ([3} where k is a fit parameter. In 
Fig.Q](b) different external fields are applied for a constant 
time. The dashed line is a guide to the eye. 



ing 40 /im saturation is obvious. 

Fig. [2] not only shows the measured strengths of the 
revealed index-of-refraction grating for different period 
lengths A, but also a calculation of the degree of modu- 
lation. Here assumptions on the domain size distribution 
have to be made: A Gaussian distribution centered at 
with a width of 5 fim is assumed. The height of the plot- 
ted curve is used as a fit parameter. The basic message is 
that if a broad domain size distribution is assumed, as it 
is found in real crystals by the PFM measurements, the 
model predicts that no revealed grating should be found 
for period lengths smaller than the length of most of the 
domains. For higher period lengths the grating quality 
increases. These predictions are in agreement with mea- 
sured data. 

To check whether a space-charge grating is revealed 
all over the crystal or only in some parts of it, the index 
grating induced by the space-charge field is imaged with 
differential interference contrast microscopy (DIC) after 
repoling the crystal. This method allows to detect very 
small gradients of the index-of-refraction in the crystals. 
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Grating period A [urn] 

FIG. 2: The index-of-refraction modulation An is plotted ver- 
sus the period length A of the fixed grating. To correct for dif- 
ferent degrees of poling, An was divided by —(1 — g) ln(l — q) 
(e.g. Eq. ©). The curve is the product of a fit parame- 
ter A and the caluclated k, which determines the degree of 
modulation Aq. A Gaussian distribution is assumed for the 
distribution of domain lengths. 




FIG. 3: A revealed grating is imaged with the differential 
interference contrast method. The period length A of the 
electrically fixed grating is 9.6 fj,m. The z-axis is aligned hor- 
izontally. 



We use this as a method to observe restored space-charge 
fields and to determine qualitatively the shape of fixed 
domain gratings. This is possible because after repol- 
ing of the crystal uncompensated charges appear at the 
former domain walls. These charges cause electro-optic 
index changes which are a replica of the original domain 
pattern [13]. Figure[3] shows the result. As it can be seen, 
the grating reappears all over the crystal. 

Piezoelectric force microspcopy (PFM) is used to di- 
rectly image the ferroelectric domain structure after elec- 
trically fixing a space-charge grating. The PFM is a 
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FIG. 4: A fixed domain grating with a period length A of 
11.7 /xm is imaged by piezoelectric force microscopy. The two 
grey values represent the two directions of the spontaneous 
polarization. The z-axis is aligned horizontally. 

scanning force microscope operated in contact mode with 
an additional alternating voltage applied to the tip [20j . 
In piezoelectric samples this voltage causes thickness 
changes and therefore vibrations of the surface which lead 
to oscillations of the cantilever that can be read out with 
a lock-in amplifier. In the present work we imaged the 
non-polar faces of LiNb03. 

The faces of the crystal are polished down stepwise 
by some hundred micrometers to ensure that volume do- 
mains are visualized. Figure[4] shows that ferroelectric 
volume domains can be found after electrical fixing of 
holograms in lithium niobate crystals. However, a peri- 
odic structure is not immediately evident. In fact, the 
length of some areas with inverted spontaneous polariza- 
tion are bigger than one period of the fixed grating. 

However, from the PFM images we conclude that elec- 
trical fixing indeed leads to ferroelectric volume domains 
with inverted spontaneous polarization in head-to-tail 
configuration, but the periodicity is not evident locally. 
The DIC images suggest that averaged over the whole 
crystal thickness a well defined grating is present. Both 
observations are consistent with the predictions of the 



model. The degree of modulation is too small to be ob- 
vious from direct domain visualization, but an average 
periodicity is visible which can be detected by the help 
of light diffraction at the revealed space charge grating. 

After revealing, local fields are created according to the 
former domain structure. Typical values of such fields 
E\ oc presumably are determined by material parameters, 
as the spontaneous polarization and field-limiting effects 
like the electric-breakdown field. The E\ oc form a corre- 
sponding index-of-refraction grating via the electro-optic 
effect which is detected by diffraction of a light beam. 
The effective, averaged field E g, which is measured by 
this method, is given by the modulation degree multi- 
plied by the local electric field -Eioc- Thus theoretical 
predictions for Aq link a material dependent factor E\ oc 
with the experimental value E c q, i.e., these values are 
proportional to each other. 

Alltogether it is found that the predictions of the model 
are in agreement with all available experimental data, in- 
dicating strong support for the model presented herein. 
None of the assumptions of the theory is specific to the 
material lithium niobate. So the stochastic description 
offers a general way to analyze and discuss experimen- 
tal dependences in all electrical fixing experiments which 
have been performed in a large variety of ferroelectric 
materials. For Example in Ref. [2l[ the same dependence 
of revealed grating strengths on the degree of poling and 
the same dependence on the grating period A is found 
for electrical fixing in barium titanate crystals. 

Furthermore, the stochastical description enables 
quantitative predictions. Provided better knowledge of 
local fields after revealing is gained, the factor j3/a, i.e. 
the influence of local fields on the poling dynamics, can 
be analyzed. This is a point which is hard to access by al- 
ternative means. As the influence of defects presumably 
can be described by local fields, deeper understanding of 
ferroelectric domain switching becomes possible. 
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Stiftung is gratefully acknowledged. 
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